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ABSTRACT 

The expression, replication and repair of eukaryotic 
genomes require the fundamental organizing unit 
of chromatin, the nucleosome, to be unwrapped 
and disassembled. We have developed a quanti- 
tative model of nucleosome dynamics which 
provides a fundamental understanding of these 
DNA processes. We calibrated this model using 
results from high precision single molecule nucleo- 
some unzipping experiments, and then tested its 
predictions for experiments in which nucleosomes 
are disassembled by the DNA mismatch recognition 
complex hMSH2-hMSH6. We found that this 
calibrated model quantitatively describes hMSH2- 
hMSH6 induced disassembly rates of nucleosomes 
with two separate DNA sequences and four distinct 
histone modification states. In addition, this model 
provides mechanistic insight into nucleosome 
disassembly by hMSH2-hMSH6 and the influence 
of histone modifications on this disassembly 
reaction. This model's precise agreement with 
current experiments suggests that it can be 
applied more generally to provide important mech- 
anistic understanding of the numerous nucleosome 
alterations that occur during DNA processing. 

INTRODUCTION 

Eukaryotic DNA is packaged into chromatin by wrapping 
around histone protein octamers to form nucleosomes 
(1,2). These nucleosomes sterically occlude replication, 
transcription and repair complexes from their DNA 
target sites. DNA processing complexes appear to access 
their target sites via nucleosome unwrapping (3-5) and 
disassembly (6). However, the mechanism(s) associated 
with these nucleosome structural alterations are not well 
understood. An essential component in appreciating these 



nucleosome alterations is a quantitative and predictive 
model describing DNA unwrapping and disassembly 
from the histone octamer. Such a model will provide a 
critical tool in understanding in vivo DNA processing 
within chromatin. 

The key concept in a model of nucleosomal DNA 
unwrapping is the unwrapping free energy landscape. 
This landscape quantitatively captures the free energy 
cost of unwrapping each individual base pair from the 
histone octamer. Previous studies determined qualitative 
features of the structure of this free energy landscape (4,5). 
Here, we develop a free energy landscape that not only 
accounts for these quahtative observations, but is 
cahbrated from precise mechanical unzipping measure- 
ments of nucleosomal DNA (7). These experiments 
apply a force to puU apart the DNA base pairing, which 
concurrently unwraps the DNA from the histone octamer. 
We extract detailed information about the energy land- 
scape from these experiments. By carefully adjusting our 
energy landscape, we obtain an excellent agreement with 
the highly complex measured unzipping profiles. 

We then use this cahbrated free energy landscape to 
model a completely separate experiment that investigates 
nucleosome disassembly by the human DNA mismatch 
recognition complex, hMSH2-hMSH6. This complex 
appears to function as a molecular switch that specifically 
binds DNA mismatches in the ADP bound form (8,9). 
Upon the exchange of ADP for ATP, hMSH2-hMSH6 
forms a sliding clamp that diffuses along duplex DNA. 
This allows for iterative loading of multiple hMSH2- 
hMSH6 complexes onto the duplex DNA surrounding 
a mismatch, which directs the excision reaction during 
DNA mismatch repair. 

This mismatch recognition process must occur within 
chromatin. Recently, it was demonstrated that hMSH2- 
hMSH6 disassembles nucleosomes near a DNA mismatch 
and that histone modifications in the DNA-histone inter- 
face of the nucleosome dramatically enhance the rate of 
disassembly (10,11). These results indicate that hMSH2- 
hMSH6 not only recognizes DNA mismatches but 
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facilitates tlie DNA excision reaction by ensuring the 
DNA surrounding a mismatch is nucleosome and 
perhaps protein free. It was proposed that hMSH2- 
hMSH6 disassembles nucleosomes by trapping DNA 
unwrapping fluctuations through iterative loading of 
hMSH2-hMSH6 at the DNA mismatch (10). 

We modeled this nucleosome disassembly process with 
our free energy landscape and found that disassembly of 
unmodified nucleosomes can be quantitatively understood 
using the same free energy landscape applied to unzipping 
experiments. We explain the influence of histone post 
translational modifications (PTMs) and DNA sequence 
on the disassembly rates with modest changes to the 
free energy landscape based on the location of the modi- 
fications and the known influence of DNA sequence and 
histone modifications on DNA-histone binding free 
energies (4,11,12). Furthermore, our model provides 
detailed mechanistic insight into the process of nucleo- 
some disassembly by hMSH2-hMSH6. 



MATERIALS AND METHODS 

Calculation of the dynamics of nucleosome unwrapping 

We treat the time evolution of the probabihty that n 
base pairs of nucleosomal DNA are unwrapped as a 
continuous-time Markov process. This allows us to write 
the rates of wrapping ktemap(n) and unwrapping /cunwrap(«) 
in a rate matrix r„_„', which is defined by 



'^n,ti+\ — ^unwrap('^) 



rewrap 



(n), 



(1) 
(2) 



with all other entries in T„^„/ equal to zero. This rate matrix 
satisfies the master equation 



dt 



(3) 



where P„{t) is the time-dependent probability of observing 
the unwrapped state n. The rate matrix T,,,,,' contains all 
of the information about the influence of external con- 
straints, such as an applied force or a protein complex, 
on the probability and time evolution of each partially 
unwrapped nucleosome state n. Below, we describe how 
we determine the rate matrix that describes the influence 
of a mechanical unzipping force and a pressure from 
bound hMSH2-hMSH6 complexes on the nucleosome 
unwrapping probability and its time evolution. However, 
once this rate matrix is set up, the following solution to 
the master equation may be apphed to any system where 
nucleosomes are unwrapped. 

We solve the eigenvalue problem for this master 
equation to obtain P„{t). We compute the eigenvalues 1,- 
and eigenvectors F,_„ of J^j^ numerically (13). We then use 
the inverse of the matrix of eigenvectors F^,' to express the 
initial state in the basis of the eigenvectors 



(4) 



where Pn{0) is the initial state in the basis of the 
unwrapping states. If the nucleosome is fully wrapped at 
time zero, then Po(0) = 1 and ^,(0) = 0 for / 7^ 0. This 
then allows us to write the solution to Equation 3 as 



^«(0 = ^^"v,A(0)exp(-|l,|0- 



(5) 



The term exp(— |A,|/) propagates each eigenvector with 
a decay rate given by its eigenvalue, and we multiply by 
F,;„ to transform back into the basis of unwrapping states. 
This solution allows us to calculate the complete dynamics 
of nucleosome unwrapping rapidly. 

Rates of nucleosome unzipping 

In unzipping experiments, an applied force F breaks 
apart the base pairing of DNA containing a nucleosome, 
which simultaneously unwraps the nucleosome. We 
compute A;unwiap(«) and k^em-Apin), the rates of nucleosome 
unwrapping and rewrapping at n unwrapped base pairs, 
using the Gibbs free energy difference between the change 
in our energy landscape Gnuc(«) due to unwrapping and 
the work done on the nucleosome by the applied force. 
We multiply the base rate for wrapping fluctuations /cq by 
a Gibbs factor, resulting in 



t („\ _ r, „[(f-fu„zip("))aAA--AG„uc(")]Asr 
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(6) 
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where kg is Boltzmann's constant and T is the tem- 
perature. The applied force F acts over a distance Ax 
per opened base pair, and Funupin) is the force required 
to break the DNA base pairing and stretch the DNA. Of 
the remaining force, a fraction a increases the unwrapping 
rate and (1 — a) decreases the rewrapping rate. The par- 
ameter a represents the relative position of the transition 
state between n and n+\ unwrapped base pairs. We may 
then insert these rates of unwrapping and rewrapping into 
Equations 1 and 2, and follow the procedure outhned 
in the previous section to calculate the dynamics of 
nucleosome unwrapping in unzipping experiments. 

Calculation of dynamics in the presence of 
hMSH2-hMSH6 

When hMSH2-hMSH6 is present we must model its binding 
to the DNA, the interactions between hMSH2-hMSH6 
molecules, and the interactions between hMSH2-hMSH6 
and the nucleosome. Given that the rate of hMSH2- 
hMSH6 diffusion along the DNA is faster than all other 
rates in the system, we treat hMSH2-hMSH6 diffusion as 
being in equilibrium with respect to all other processes in the 
system. We may then model the hMSH2-hMSH6 molecules 
as a ID gas of hard spheres. This is known as a Tonks gas, 
and has a partition function given by (14) 



ZNj = {l-NGf/N\, 



(8) 



where N is the number of hMSH2-hMSH6 molecules 
bound to the DNA, ct is the footprint of hMSH2-hMSH6 
in base pairs and / is the length of the DNA segment avail- 
able to hMSH2-hMSH6 in base pairs (Figure 2A). 
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To confirm that use of the Tonks gas solution is vahd, for 
select cases we also performed simulations of the dynamics 
including each hMSH2-hMSH6 particle individually. We 
find that approximating the hMSH2-hMSH6 behavior as 
a Tonks gas gives results nearly identical to these simula- 
tions, and is much more computationally efficient 
(Supplementary Figure SI). 

We can then use the Tonks gas partition function to 
calculate the probability that the bound hMSH2-hMSH6 
molecules will block nucleosome rewrapping or the 
binding of more hMSH2-hMSH6. The probability that 
there are no hMSH2-hMSH6 particles within x base 
pairs of one end of the available DNA segment for 
X < (/ - Ncj) is 



l-Na-. 
1-Ng 



(9) 



To calculate the rates of unwrapping, A:unwrap(«)5 and 
rewrapping, ^iewrap(n), we multiply the base rate for 
wrapping fluctuations ko by a Gibbs factor determined 
from our free energy landscape, as in Equations 6 and 7 



k 



unwrap 



(«) = ^-oe-^G„„,(„)/i.-^r 



kiemap(n) = ^'0Pa',/(1), 



(10) 

(11) 



where Pnj(I) is the probabihty that an hMSH2-MSH6 
particle does not block rewrapping of one base pair. 
Similarly, the rate of hMSH2-hMSH6 binding at the 
mismatch is 



^^nT ~ ^'(),?''pa',/(«m+o'/2), 



r bind 



(12) 



where /cq'"'* is the rate of binding to bare DNA of length 
/, riM is the mismatch position, and pj^ /{hm+ct/T) is the 
probability that binding is not blocked by the presence 
of another hMSH2-hMSH6 molecule. This assumes that 
the mismatch is within 1.5ct base pairs of the end of the 
DNA segment, such that there cannot be any hMSH2- 
hMSH6 molecules between the end of the DNA and the 
mismatch when an hMSH2-hMSH6 molecule is bound 
at the mismatch. If this is not true, one must sum over 
ah combinations for the number of hMSH2-hMSH6 
complexes bound on either side of the mismatch. 

We then calculate the dynamics of nucleosome unwrap- 
ping foUowing the method described in the 'Calculation of 
the dynamics of nucleosome unwrapping' section. Instead 
of indexing the probabihty in Equation 3 based only on «, 
we use II = NL+n, where is the number hMSH2- 
hMSH6 molecules already bound to the DNA, and L is 
the maximum number of bases which can unwrap from 
the histone before it dissociates from the DNA. This 
allows us to write a rate matrix Tj,,,,/ which includes not 
only rates of nucleosome wrapping and unwrapping, but 
also hMSH2-hMSH6 binding and dissociation. 

RESULTS 

Mechanical unzipping of nucleosomes 

The most detailed method yet devised to probe nucleo- 
some dynamics is mechanically unzipping the base 



pairing of DNA containing a nucleosome (7). In experi- 
ments performed by the Wang group (7), the nucleosome 
is positioned on the 601 high-affinity sequence (15). One 
strand of the DNA molecule is attached to a surface, and 
the other strand to a bead confined in an optical trap 
(Figure lA). By carefully controlling the movement 
of the optical trap, a constant force may be exerted 
on the DNA. This force sequentially breaks the 
DNA base pairing, increasing the contour length of the 
single-stranded DNA between the bead and the surface. 
Measuring the position of the bead as a function of time 
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Figure 1. A free energy landscape consistent with measurements of 
mechanical unzipping of nucleosomes. (A) Experimental setup for 
unzipping experiments. One strand of the DNA is attached to a 
surface, and a constant force F of 28 pN is applied to the other end 
via a bead in an optical trap, slowly pulling apart the DNA base 
pairing. As this unzipping of the base pairing moves through the nu- 
cleosome, the DNA must also be unwrapped from the nucleosome. 
(B) Dwell time as a function of the number of DNA base pairs (bp) 
unwrapped from the nucleosome. The solid line is determined from 
mechanical unzipping experiments, while the dashed line is calculated 
from (C) the free energy landscape for DNA unwrapping from a 
nucleosome on the 601 positioning sequence. Significant pauses when 
20-30 bp are unwrapped (within the entry-exit region) and when 
60-80 bp are unwrapped (within the dyad region) correspond to areas 
where the free energy landscape changes rapidly. 
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establishes the number of DNA bases that are unpaired. 
When DNA wrapped around the nucleosome becomes 
critically unpaired, it must dissociate from the nucleo- 
some. Thus, the number of broken base pairs equals the 
number of base pairs unwrapped from the nucleosome 
(sohd line. Figure IB). This method probes the interaction 
between the nucleosome and the DNA as a function of the 
number of base pairs unwrapped from the nucleosome. 

We calculate an energy landscape for unwrapping 601 
positioning sequence DNA from a nucleosome, which is 
consistent with this unzipping data (Figure IC). To obtain 
this landscape, we first create a trial landscape that we use 
to compute rates of nucleosome unwrapping. We then use 
these rates to calculate unzipping dynamics, and compare 
the result with experiments. Finally, we iterate the trial 
landscape until we obtain good agreement with experi- 
ments (Figure IC) while respecting several quahtative 
features of the landscape known before our study. 

In order to calculate experimental observables given an 
energy landscape, we compute A'unwrap(«) and A'rewiap(«), 
the rates of nucleosome unwrapping and rewrapping 
at 77 unwrapped base pairs, using Equations 6 and 7. 
In these experiments, the total applied force F is held 
constant at 28 pN (7). The force Funzip(«) needed 
to unzip the base pairing and stretch the DNA can be 
calculated as a function of the DNA sequence 
(Supplementary Figure S2) (16). This results in an effect- 
ive force responsible for nucleosome unwrapping 
averaging 10.9 pN, and ranging from 8.6 pN to 12.8 pN. 
Similarly, the increase in contour length from opening one 
base pair A.y can be calculated to be 1. 12 nm at this force 
(16). Since the relative position a of the transition state 
between n and n+l unwrapped base pairs is not known 
under the geometry of the experiments, we treat it as a fit 
parameter. 

We find that our results are not very sensitive to the 
value of ko, the rate of nucleosome rewrapping in the 
absence of an applied force. We can fit the unzipping 
data with values of between 10^ to 10^ base pairs per 
minute by simultaneously making small adjustments 
to the nucleosome energy landscape and a. We used the 
value of A'o = 6 X 10^ base pairs per minute that is in the 
middle of this range for all subsequent calculations. 

The trial landscape must conform to a number of 
experimental observations. First, the free energy land- 
scape should be steepest in the dyad region (between 
about 63-84 bp), since the nucleosome contains the 
strongest contacts in this region (2,4). Second, the total 
height of the free energy landscape must be no more 
than 42 keT, of the free energy cost to disassemble a 
nucleosome from a compacted chromatin fiber in 
Xenopus egg extracts (17). The buffer conditions and 
presence of histone HI in Xenopus egg extracts allow 
for strong interactions between nucleosomes, and promote 
chromatin compaction. We expect the free energy cost 
to disassemble the recombinant histones examined here 
to be significantly smaller, as the measurements are per- 
formed under conditions where chromatin compaction 
does not occur. 

Third, the steps in the landscape occur approximately 
every 5.25 bp, which is the average spacing between DNA 



contacts in the nucleosome (2) and between peaks in 
the unzipping measurements (7). The positions of the 
peaks in the unzipping data (7) are extremely sensitive 
to the spacing between the steps in the unzipping land- 
scape. For example, if we force all steps to be separated 
by exactly 5.25 bp, we are not able to accurately fit 
the unzipping data (Supplementary Figure S3). We can 
thus directly read off the spacing of the steps from the 
positions of the peaks in the unzipping data (7). This 
results in shght variations of the step spacing around 
the average of 5.25 bp consistent with the nucleosome 
crystal structure, which shows that DNA contacts in 
the nucleosome are not evenly spaced (2) either. We use 
discrete steps for simplicity, as we find no signifi- 
cant difference from using more gradual steps. Lastly, 
the landscape cuts off shortly after the dyad region, 
such that the histone octamer irreversibly dissociates 
from the DNA and cannot bind again once the system 
reaches this cut off This corresponds with the observa- 
tion that nucleosomes disassemble after the dyad is 
disrupted (7,18). 

The free energy landscape is varied until the calculated 
DNA unzipping dynamics (dashed line. Figure IB) 
agree with the experimental data. The free energy land- 
scape shown (Figure IC) conforms to all the constraints 
outlined above and also accurately reproduces the 
measured unzipping dynamics. Thus, we are able to 
extract detailed information about the structure of the 
nucleosome free energy landscape from these studies. 

Because we do not know the value of a, the landscape is 
not completely constrained by the unzipping experiments. 
By making relatively small adjustments to the value of a, 
we can obtain free energy landscapes that fit the unzipping 
data about as well as the one shown (Figure IC), but that 
have total height as small as 12 kgT oi as large as 27 kjjT. 
We select a landscape with total height of 23.8 kgT as 
this allows us to simultaneously fit measurements of 
nucleosome dissasembly by hMSH2-hMSH6 described in 
the next section. 

Nucleosome disassembly by hMSH2-hMSH6 

To demonstrate the utility of our nucleosome free energy 
landscape, we apply it to measurements of the rate of 
nucleosome disassembly induced by the DNA mismatch 
recognition complex, hMSH2-hMSH6 (10,11). In these 
experiments, one end of a duplex DNA molecule is 
labeled with biotin, while the other end contains a nucleo- 
some positioning sequence wrapped into a nucleosome 
(Figure 2A). The DNA molecule contains 99 or 75 base 
pairs between the biotin and the 5S or 601 positioning 
sequence, respectively, and a mismatch 21 base pairs 
from the biotin labeled end. hMSH2-hMSH6 binds the 
DNA mismatch and following ATP binding forms a 
clamp that slides freely along the DNA (9,19). 
Streptavidin bound to the biotin labeled end of the 
DNA prevents hMSH2-hMSH6 from sliding off that 
DNA end. At the opposite end, hMSH2-hMSH6 disas- 
sembles the nucleosome, allowing for hMSH2-hMSH6 
to slide off the DNA. 
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Figure 2. Model predictions compared with experimental meas- 
urements for nucleosome displacement by hMSH2-hMSH6. (A) 
Graphical depiction of the experimental system [10]. Streptavidin 
bound to the biotinylated left end of the DNA prevents hMSH2- 
hMSH6 (MSH) from sliding off this end, while a nucleosome inhibits 
hMSH2-hMSH6 from sliding off the right end of the DNA. There is 
a mismatch located «m = 21 bp from the biotin bead, at which 
hMSH2-hMSH6 may bind. The total length of the DNA segment is 
L. The hMSH2-hMSH6 is confined to length /, the DNA between the 
biotin bead and the nucleosome. The footprint of one hMSH2-hMSH6 
molecule is a = 25 bp. Free energy landscapes for unwrapping nucleo- 
somes on (B) the 601 positioning sequence and (C) the 5S positioning 



Nucleosome disassembly by hMSH2-hMSH6 was 
demonstrated to depend on ATP binding but not ATP 
hydrolysis (10). We therefore consider a model where no 
ATP is hydrolyzed and hMSH2-hMSH6 cannot directly 
exert a force on the nucleosome. Rather, hMSH2-hMSH6 
takes advantage of fluctuations in nucleosonial DNA 
unwrapping (3,5,20), which we calculate from our 
cahbrated free energy landscape. hMSH2-hMSH6 binds 
at the DNA mismatch, forms a sliding clamp and then 
freely diffuses along the DNA. An additional hMSH2- 
hMSH6 complex can bind at the mismatch and form 
a sliding clamp once hMSH2-hMSH6 diffuses off the 
mismatch. As muhiple hMSH2-hMSH6 clamps load 
onto the DNA, they can occupy unwrapped nucleosonial 
DNA, which prevents rewrapping thereby increasing the 
fraction of unwrapped nucleosomal DNA. 

To accurately model hMSH2-hMSH6, we use known 
rates of hMSH2-hMSH6 shding and binding. A rate of 
5 per minute was used for hMSH2-hMSH6 binding and 
sliding clamp formation at a G/T mismatch, while 0.1 per 
minute was used for the rate of dissociation of the sliding 
clamp form of hMSH2-hMSH6 from DNA (8,19,21,22). 
A rate of 3.7 x 10^ base pairs per minute was used for the 
step rate of hMSH2-hMSH6 clamps. This corresponds to 
the reported lower limit for the hMSH2-hMSH6 clamp 
diffusion constant on duplex DNA of 0.036 (im^/s (21). 

We calculated rates of hMSH2-hMSH6 induced nucleo- 
soiTie disassembly for comparison to rates measured by 
electrophoretic mobihty shift assays (10,11). For compu- 
tational efficiency, we used the fact that hMSH2-hMSH6 
stepping along the DNA is the fastest process in our 
model, which allows us to model hMSH2-hMSH6 
clamps as an equilibrated gas (14). This gas consists of 
particles of size a = 25 nm (8), which diffuse freely in 
one dimension along the DNA, but cannot pass through 
one another (Figure 2A) (see 'Materials and Methods' 
section). We find the calculated nucleosome disassembly 
rates using this approximation to be nearly identical to 
those from a full simulation including every hMSH2- 
hMSH6 complex explicitly (Supplementary Figure SI). 

DNA sequence dependence of nucleosome disassembly 
by hMSH2-hMSH6 

We begin by examining nucleosome disassembly measure- 
ments with the 601 positioning sequence. This allows us to 
employ the same free energy landscape used to model 
unzipping experiments (solid hne. Figure 2B). We calcu- 
late the characteristic time for nucleosome disassembly 
to be 970 min, which agrees with the experimental obser- 
vation that hMSH2-hMSH6 does not measurably dis- 
assemble nucleosomes in 60 min (diamonds. Figure 2D) 
(11,23). 



Figure 2. Continued 

sequence. Numerical values for these energy landscapes are provided as 
a supplemental spreadsheet. A comparison of the results of our 
modeUng (lines) using these energy landscapes to experimental meas- 
urements (points) of nucleosome displacement by liMSH2-hMSH6 
shows excellent agreement for both (D) the 601 positioning sequence 
and (E) the Xenopus 5S positioning sequence. 
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Next, we consider the influence of the 5S nucleosome 
positioning sequence on nucleosome disassembly by 
hMSH2-hMSH6. We assume that the change in DNA 
sequence only influences the free energy landscape of nu- 
cleosome unwrapping. Therefore, we kept all model par- 
ameters fixed except for the nucleosome unwrapping free 
energy landscape. While detailed DNA unzipping meas- 
urements have not been reported for a 5S sequence, we 
expect this landscape to have the same general form as the 
601 landscape, but a smaller total height. The sea urchin 
5S sequence lowers the free energy of nucleosome forma- 
tion by 4.7 kgT relative to 601 (15) and thus we expect a 
reduction of several kgT for the Xenopus 5S sequence as 
well (10). Furthermore, the DNA sequence near the dyad 
region appears to be responsible for the enhanced DNA- 
histone binding of the 601 sequence (4,15). 

We modified the 601 landscape by decreasing it by 
a total of 2.2 k^T in the dyad region (sohd fine. 
Figure 2C). Nucleosome disassembly by hMSH2- 
hMSH6 is not sensitve to the detailed structure of the 
free energy landscape, so our choice to modify the land- 
scape at the dyad region is based on empirical evidence 
that the increased binding affinity of the 601 sequence 
relative to 5S is due mainly to contacts in the dyad 
region (4,15). This landscape gives a good fit to the experi- 
mental data for unmodified nucleosomes containing 
the Xenopus 5S sequence (10) (diamonds. Figure 2E). 
This landscape allowed us to explain the behavior of nu- 
cleosome disassembly in the absence of a mismatch 
(squares. Figure 2E) where hMSH2-hMSH6 is permitted 
to bind at any position within the DNA at a rate of 
9 X 10""* per minute, the experimental rate of hMSH2- 
hMSH6 binding to duplex DNA without a mismatch 
(8,19,22). These results suggest that modest changes 
in the free energy landscape due to alterations in the 
nucleosomal DNA sequence can lead to significant 
changes in the hMSH2-hMSH6 induced nucleosome 
disassembly rate. 

The influence of histone post-translational modiflcations 
of nucleosome disassembly by hMSH2-hMSH6 

Histone PTMs within the DNA-histone interface reduce 
the free energy of nucleosome formation (1 1,12,24,25) and 
facihtate nucleosome disassembly by hMSH2-hMSH6 
(10,11). We incorporated the influence of histone PTMs 
into our model by assuming that they specifically influence 
the DNA unwrapping free energy landscape close to the 
location of the PTM. 

We first investigated the influence of histone PTMs near 
the nucleosome dyad: the phosphorylation of histone H3 
at Til 8 [H3(T118ph)] and the double acetylation of H3 at 
K115 and K122 [H3(K115ac,K122ac)]. These modifica- 
tions occur in vivo (26) and influence DNA repair (27). 
Recently, it was reported that H3(T118ph) reduces the 
DNA-histone binding free energy by 3.5±0.3 ksT and 
that it increases the rate of hMSH2-hMSH6 induced nu- 
cleosome disassembly by 25 times (11). We found that 
a 3.8 hgT reduction of our 601 free energy landscape 
in the dyad region results in excellent agreement between 



the model predictions and the measured data (squares. 
Figure 2D). 

The dyad acetylations, H3(K1 15ac,K122ac), reduce the 
DNA-histone binding free energy by 0.8 ± 0.4 hgT (12) 
and increase the rate of nucleosome disassembly by 
hMSH2-hMSH6 (10). We found that decreasing the free 
energy landscape for unmodified 5S by 1.2 ksTin the dyad 
region fits the experimental measurements for these modi- 
fications (circles. Figure 2E). 

Histone acetylation in the DNA entry-exit region of the 
nucleosome, H3(K56ac), reduces the free energy of the 
nucleosome binding to DNA (24) and increases DNA 
unwrapping (25). From this increase in unwrapping, we 
estimate that H3(K56ac) reduces the binding free energy 
by 0.5 ±0.3 ksT (25). This modification occurs within 
newly replicated DNA and appears important for DNA 
repair (6). Furthermore, the acetyllysine mimic at this 
location, H3(K56Q), enhances the rate of nucleosome dis- 
assembly by hMSH2-hMSH6 (10). We incorporated 
H3(K56Q) into our model by reducing the energy land- 
scape by 0.5 ksT in the entry-exit region compared to the 
unmodified 5S landscape (dashed line. Figure 2C). This 
results in an excellent fit to the hMSH2-hMSH6 displace- 
ment data for H3(K56Q) (triangles. Figure 2E). 

Finally, using the energy landscapes described above, 
we reproduce the disassembly rates without a mismatch 
for H3(T118ph) (inverted triangles. Figure 2D) and 
H3(K1 15ac,K122ac) 5S nucleosomes (inverted triangles. 
Figure 2E). We use the same hMSH2-hMSH6 binding 
conditions used to describe unmodified nucleosomes 
with no mismatch on 601, i.e. without fitting any other 
parameters. As before, measurements without mismatches 
were modeled by allowing hMSH2-hMSH6 to bind at any 
position within the DNA. 



DISCUSSION 

To our knowledge, we propose the first detailed quantita- 
tive free energy landscape for unwrapping DNA from a 
histone octamer. We develop this landscape from high 
resolution nucleosomal DNA unzipping experiments (7). 
We then demonstrate that this landscape accounts for the 
different hMSH2-hMSH6 induced nucleosome disassem- 
bly rates for two separate DNA sequences and four 
separate histone modification states. All parameters are 
held constant except for the free energy landscape. To 
model modifications and DNA sequence differences we 
only introduce modest landscape alterations located 
in the regions of the landscape associated with the 
modification or sequence difference. The magnitudes of 
these alterations are consistent (Table 1) with changes 
in DNA-histone binding free energies known from experi- 
ments (11,12,15,25). These results demonstrate that the 
nucleosome unwrapping free energy landscape provides 
an important tool for understanding enzyme provoked 
nucleosome disassembly. 

The model's success in explaining numerous experimen- 
tal results indicates that it captures the essential compo- 
nents of nucleosome disassembly by hMSH2-hMSH6 
and that we can use the model to provide mechanistic 
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Table 1. Binding energy for modified nucleosomes relative to 
wild-type nucleosomes 



Sequence 


Modification 




A F ^ 
i-^J-- landscape 






(ksT) 


(ksT) 


601 


H3(T118pli) 


3.5 ±0.3 


3.8 


5S 


H3(K56Q) 


0.5 ±0.3 


0.5 


5S 


H3(K115ac,K122ac) 


0.8 ±0.4 


1.2 



"Experimentally measured binding energies [11,12,25]. 
''Binding energies determined from our free energy landscapes. 



insight into this disassembly process. The observation that 
multiple hMSH2-hMSH6 clamps loaded onto DNA 
behave as a ID gas, also known as a Tonks gas (14), 
provides a mechanistic picture of hMSH2-hMSH6 
induced nucleosome disassembly. The hMSH2-hMSH6 
Tonks gas exerts a ID pressure on the nucleosome 
(Supplementary Figure S4). The pressure builds as add- 
itional hMSH2-hMSH6 complexes load onto the DNA, 
which shifts the equihbrium amount of DNA wrapped 
into the nucleosome toward more unwrapped states. 
This enables the nucleosome to eventually reach an 
unwrapped state that allows for histone octamer 
disassociation. 

The number of hMSH2-hMSH6 clamps required to dis- 
assemble the nucleosome is three to four (Figure 3 and 
Supplementary Figure S5). However, this number of 
clamps is not sufficient to trap DNA unwrapping states 
out to the dyad. Instead, the plateau in the free energy 
landscape around 50 bp allows clamps to trap states with 
about 50 bp of unwrapped DNA often enough that an 
unwrapping fluctuation to the dyad occurs before the 
outermost clamp shdes away from the nucleosome 
(Supplementary Figure S6). 

The model indicates that the mechanism by which dyad 
modifications, H3(T118ph) and H3(K115ac,K122ac) 
influence nucleosome disassembly is distinct from the 
DNA entry-exit modification mimic, H3(K56Q). None 
of the modifications influence the distance the outermost 
hMSH2-hMSH6 clamp invades into the nucleosome. 
The dyad modifications increase the probability for an 
unwrapping fluctuation that allows for histone octamer 
release when the outermost hMSH2-hMSh6 clamp 
has trapped about 50 bp of DNA (Supplementary 
Figure S6). In contrast, H3(K56Q) increases the probabil- 
ity that hMSH2-hMSH6 has trapped an unwrapping fluc- 
tuation of about 50 bp. This increases the time period the 
nucleosome exists in a partially unwrapped state and in 
turn the probabihty that the nucleosome will unwrap to 
the dyad and release the histone octamer. These mechan- 
istic insights elucidate how fractional changes in the nu- 
cleosome unwrapping free energy landscape result in order 
of magnitude changes in the rate of disassembly. 

The time for nucleosome disassembly varies from 
about 40min to 16 h. Yet, all the processes involved 
with disassembly, including nucleosome unwrapping, 
hMSH2-hMSH6 shding, and hMSH2-hMSH6 binding, 
are significantly faster. Our model explains this apparent 
discrepancy. The first and second hMSH2-hMSH6 clamps 
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Figure 3. Distribution of the number of bound hMSH2-liMSH6 com- 
plexes needed to displace nucleosomes from the 601 positioning 
sequence. The solid line shows the result for unmodified nucleosomes, 
and the dashed fine for the dyad modification H3(T118ph). 



load on the scale of minutes. However, to load a third 
clamp, both loaded hMSH2-hMSH6 clamps must simul- 
taneously move toward the nucleosome to provide enough 
space for the third hMSH2-hMSH6 clamp. This is 
strongly disfavored by entropy. Loading even larger 
numbers of hMSH2-hMSH6 clamps requires the coordin- 
ation of nucleosome unwrapping with simultaneous 
shding of multiple clamps, which is further disfavored by 
entropy. As a consequence, the number of hMSH2- 
hMSH6 clamps required for disassembly is three to four 
and the rates for loading the fourth clamp are similar to 
the rates for nucleosome displacement (Supplementary 
Figure S7). We conclude that the binding of multiple 
hMSH2-hMSH6 clamps is the rate-limiting step in 
nucleosome displacement. 

Our model not only quantitatively explains recent ex- 
periments by Javaid et al. (10) and North et al. (11), but 
also explains the apparently contradictory studies that 
report nucleosomes are a barrier to hMSH2-hMSH6 
shding (23,28). The studies by Li et al. (23) observed no 
nucleosome disassembly by hMSH2-hMSH6. While their 
construct was similar to studies by Javaid et al. and North 
et al., Li et al. only used the 601 positioning sequence 
and unmodified histone octamers. As discussed above, 
we find that the 601 sequence largely prevents disassembly 
by hMSH2-hMSH6 with unmodified histone octamer, 
which was confirmed by North et al. (11). Gorman et al. 
(28) reported that nucleosomes create a boundary to yeast 
MSH2-MSH6 diffusion. Their experiments involve single 
molecule tracking of quantum dot labeled MSH2-MSH6 
clamps on lambda DNA that contains nucleosomes but no 
mismatch. They track MSH2-MSH6 along the DNA after 
washing out free MSH2-MSH6. Because there is no 
mismatch and no additional MSH2-MSH6 to load onto 
the DNA, multiple MSH2-MSH6 clamps cannot be 
loaded onto the DNA between two nucleosomes. Thus, 
under the experimental conditions used by Gorman 
et al., the pressure of the ID gas cannot reach the values 
required to disassemble a nucleosome and the nucleosome 
becomes a barrier to the sliding of a single MSH2-MSH6 
clamp as they report. 

The successful application of our nucleosoine 
unwrapping free energy landscape indicates that it will 
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prove useful for modeling a wide variety of experiments 
involving nucleosome unwrapping, and many different 
interactions between DNA binding proteins and nucleo- 
somes in vivo. In particular, when a protein does not 
interact specifically with either the DNA or the nucleo- 
some, as is the case with hMSH2-hMSH6, our free 
energy landscape completely describes the dynamics of 
the system. Even when this is not the case, the landscape 
should prove an invaluable tool. The model will facilitate 
the separation of the DNA-histone binding from add- 
itional interactions. This in turn will allow for determin- 
ation of interactions between the nucleosome and other 
protein complexes, such as chromatin remodelers that dis- 
assemble, reposition and unwrap nucleosomes. 

SUPPLEMENTARY DATA 

Supplementary Data are available at NAR Online. 
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